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Abstract 


One version of the multichannel theory for electron-target scattering based on the 
Schwinger variational principle, the SMC method, requires the introduction of a pro- 
jection parameter (K. Takatsuka and V. McKoy, Phys. Rev., A24, 2473 (1981)). The 
role of the projection parameter a is investigated and it is shown that the principal- 
value operator in the SMC equation is Hermitian regardless of the value of a as long 
as it is real and nonzero. In a basis that is properly orthonormalizable, the matrix 
representation of this operator is also Hermitian. The use of such basis is consistent 
with the Schwinger variational principle because the Lippmann-Schwinger equation 
automatically builds in the correct boundary conditions. Otherwise, an auxiliary con- 
dition needs to be introduced, and Takatsuka and McKoy’s original value of a is one 
of the three possible ways to achieve Hermiticity. In all cases but one, a can be un- 
coupled from the Hermiticity condition and becomes a free parameter. An equation 
for a based on the variational stability of the scattering amplitude is derived; its solu- 
tion has an interesting property that the scattering amplitude from a converged SMC 
calculation is independent of the choice of a even though the SMC operator itself is 
a-dependent. This property provides a sensitive test of the convergence of the calcu- 
lation. For a static-exchange calculation, the convergence requirement only depends 
on the completeness of the one-electron basis, but for a general multichannel case, the 
a-invariance in the scattering amplitude requires both the one-electron basis and the 
N+l-electron basis to be complete. The role of a in the SMC equation and the con- 
vergence property are illustrated using two examples: e-CO elastic scattering in the 
static-exchange approximation, and a two-state treatment of the e-H 2 X 1 !]^ — + b 3 E+ 
excitation. 
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I. Introduction 


Most of the recent progress in electron-atom and electron-molecule collisions are 
based on variational principles. The Kohn 1-4 (particularly the S-matrix Kohn 5,6 ) vari- 
ational principle, the variational R-matrix method, 4,7 and the Schwinger variational 
principle 8-10 have ail met with considerable success. Among the methods based on 
the Schwinger variational principle, the multichannel formulation of Takatsuka ancL 
McKoy, 11,12 sometimes referred to as the SMC method, is unique in the sense that 
its dynamical equation is derived from a combination of Lippmann-Schwinger and 
Schrodinger equations. In this method, the partitioning of open and closed channel 
space is achieved by a projection operator which operates on the target electrons alone. 
The open channel contribution is described by a projected Lippmann-Schwinger equa~ 
tion and the closed channel contribution by the Schrodinger equation. Takatsuka and 
McKoy introduced a paramater ‘a 5 to the projection operator. They argued that the 
choice of a is uniquely determined by the requirement that the the principal-value 
operator in their equation, represented in a basis which includes non 6-function or- 
thonormalizable continuum functions, must be Hermitian. Numerical calculations 13-20 
based on this choice of a are generally in good agreement with experiment and with 
calculations using other methods. 

Recently we initiated a study of the parameter a, which we call the ‘projection 
parameter 5 , in the SMC method. We begin with the premise that an optimal choice 
of a should provide a stable cross section and then formulate an equation for a based 
on the variational stability of the scattering amplitude. Analysis of the solution of this 
equation shows that the scattering amplitude deduced from a fully converged solution 
of the SMC equation is independent of a even though the SMC operator itself is a- 
dependent. This seemingly paradoxical result arises from the fact that the role of 
a in the SMC equation is related to how the Lippmann-Schwinger and Schrodinger 
equations are combined. A fully converged solution should satisfy both the Lippmann- 
Schwinger and Schrodinger equations such that the result is independent of how the 
SMC operator is partitioned. In this respect, the variational stability of the scattering 
amplitude with respect to a furnishes a stringent test for the convergence of an SMC 
calculation. We note that in most variational calculations, such as the determination 
of scattering wave functions based on the stability of the scattering amplitude or phase 
shift, and bound state wave functions based on the stability of the total energy, the 
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solutions are generally unique. The solution of the a-equation is unusual in this respect 
because, for a converged calculation, it provides an infinite (but denumerable) number 
of solutions. 

While the study of a for fully converged calculations is a purely formal exercise, 
for calculations using a finite basis, the present result is very useful in determining 
an optimal basis and testing convergence. In this respect, the test of a-invariance 
provides a unique window to the quality of the SMC calculations that apparently has 
no counterpart in other methods. It also furnishes a ‘litmus’ test for pseudo-resonances 
because pseudo-resonances are sensitive to the value of a. Finally, for studies using 
small to medium size basis sets, the global stability of the scattering cross section with 
respect to a may serve to expedite the convergence of the calculation. 

In Sec. II, we review the SMC method and the original choice of a by Takat- 
suka and McKoy 11 ’ 12 based on the Hermiticity requirement of the principal- value SMC 
operator. We show that the principal-value SMC operator is Hermitian for arbitrary, 
real, and nonzero a. The matrix representation of the principal-value SMC operator 
is also Hermitian provided that the basis functions consist of either L 2 functions and 
<5-function orthonormalizable continuum functions, or purely L 2 functions. However, 
if the continuum functions used are not ^-function orthonormalizable, Takatsuka and 
McKoy’s choice of a is one way to ensure Hermiticity. Hermiticity can also be achieved 
by two other means: iterating the projected Lippmann-Schwinger equation or enforcing 
a matrix element of E — H to be zero. In Sec. Ill we discuss the difference between 
an iY-electron and an N + 1-electron projection operator and how the introduction of 
a is related to this difference. Based on the variational principle, an equation for a is 
derived in Sec. IV and the nature of its solutions is discussed. In two specific examples, 
elastic scattering in the static-exchange approximation and multichannel close-coupling 
calculations, the a-invariance is demonstrated by substituting converged wave functions 
from other methods into the SMC equation. Also, for the case of the lowest pole of H“, 
we show that the SMC equation reduces to the Schrodinger equation and the solution 
is again independent of a. We then propose a search procedure to locate the region of 
a-invariance and establish a stable solution of the scattering amplitude. Two numerical 
examples are presented in Sec. V to illustrate the nature of the a-variation and the 
type of insight it can provide. Sec. VI concludes the paper. 
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II. The projection parameter in the SMC method 
The Hamiltonian for the electron + target system is given bv 


H — II n + T/v+i + V. (1) 

where H /v is the target Hamiltonian, XW+] is the kinetic energy operator for the con- 
tinuum electron and V the interaction potential. In the SMC formulation, a projection 
operator P, originally introduced by Takatsuka and McKoy, 11 defines the open-channel 


space in terms of the eigenfunctions <f> m 

of the target Hamiltonian, 


M 



p=Y, 1 * 41 , 2 ,. 

m= 1 

W)|, 

(2) 

and 



H N \$ m 

) = 

(3) 


The summation in Eq. (2) is over all energetically accessible target states. Notice that P 
is an iV-electron projection operator, unlike the Feshbach projection operator 28 which 
is an N -f 1-electron projection operator. In using P to partition the N+ 1-electron 
space, Takatsuka and McKoy multiply P by a projection parameter a. 

*<+> = aP*(+> + (1 — aP)¥+\ (4) 

Here is an iV-f- 1-electron wave function. At this point the parameter a is arbitrary. 
Its variational determination is the subject of the present investigation. The dynamical 
equation in the ai^-space is obtained from the projected Lippmann-Schwinger equation, 

aP¥ n +) = aPS n + a(? ( p +) Vtf <+>, ( 5 ) 

where S n is the eigenfunction of the non-interacting Hamiltonian, H ^ + T^+i- The 
projected Green s function, (?p ^ , with outgoing spherical wave boundary conditions, 
is defined in the open channel space by, 


G ( p +) 


1 V' 1* \ ex Pi ik m\rN+i - r N+1 'I) 

?.7T 2^ I m ' ir. r„ , . M 


rN+1 — rw+1 '\ 


( 6 ) 


The contribution from the complementary space is deduced from the Schrodinger equa- 
tion. 


(1 -aP)H¥ n +) = 0. 


(7) 
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H = E-H. 


(8) 


The SMC equation is derived 11 by combining Eqs. (5) and (7) and multiplying by 1/a, 

A (+) « ( n +) = VPS n , (9) 

where the SMC operator A^ with outgoing spherical wave boundary condiions is given 

by 

Al+> = ( hpv + VP) - VG { + ] V + -H - \{PH + (10) 

I A a Z J 

Lima and McKoy 29 have demonstrated that the solution of Eq. (9) provides the com- 
plete wave function. 

The Schwinger functional for the scattering amplitude, constructed from Eq. (9) 
is 

fmn = -£Wl*‘ + >) + {*l“ ) |V'|5n> - l^ ( + ) |* i n +1 ) j . (11) 

The principai-value SMC operator, A, is obtained by replacing in Eq. (10) 

with the projected principal-value Green’s function. Similarly, the wave functions in 
Eqs. (7), (9), and (11) should be replaced by those using standing wave boundary con- 
ditions. The principal-value operator was used in the original derivation of Takatsuka 
and McKoy, but most subsequent numerical calculations used 

Takatsuka and McKoy, 11,12 and Lima and McKoy 29 argued that the variational 
stability of the Schwinger functional / mn requires the principal- value SMC operator A 
to be Hermitian, or equivalently, A^^ = A^“\ While this condition is readily satisfied 
for trial functions consisting only of L 2 functions, Takatsuka and McKoy noted that 
if the trial functions included (shielded) spherical Bessel and Neumann functions, the 
representation of A in this basis is generally non-Hermitian because the Hamiltonian 
matrix in non-Hermitian. In this case, the Hermiticity of the A-matrix is retrieved if 
the projection parameter is chosen to be 

a = N + l. (12) 

Equation (12) guarantees the following relationship to hold for the continuum portion 
of the open-channel functions, 

{PV m \-H - hpH + HP)\PV n ) = 0. (13) 

a Z 
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Since the non-Hermitian part of the A-matrix involving continuum, functions is iden- 
tically zero, Eq. (12) ensures the Hermiticity of A in a basis which includes spherical 
Bessel and Neumann (or Hankel) functions. 

Two aspects of Eq. (12) deserve further consideration. We note that Eq. (13) is 
also satisfied if 

<PT m |f/|PT n ) = 0. (14) 

Thus even if the continuum functions in the basis results in a non-Hermitian iT-matrix, 
the Hermiticity of A can be achieved using conditions other than Eq. (12). This point 
will be further considered in Sec. IV. 

Another aspect of Takatsuka and McKoy’s result is that, if the Hermiticity of A 
is used to determine a, the choice of a is basis set dependent and Eq. (12) is nec- 
essary if the basis set employed gives a non-Hermitian ff-matrix. We note that the 
Hamiltonian operator , corresponding to a physical observable, must be Hermitian. 30 
Furthermore, the Hermiticity of an operator depends only on the fact that it represents 
an observable, but not on its eigenvalue being in the discrete or continuum part of the 
spectrum. However, the continuum eigenfunction of a Hermitian operator must be 
5-function orthonormalizable. 31 In order to retain the Hermiticity in the matrix rep- 
resentation of H , the basis set used needs to be orthonormalizable either in terms of 
Kronecker 5’s (L 2 functions) or Dirac 5-functions (continuum functions). However, in 
scattering calculations using Kohn or R-matrix methods, the phase shifts or T-matrices 
are deduced from the wavefunction at the asymptotic boundary. It is necessary to in- 
clude both (shielded) spherical Bessel and Neumann functions in the trial basis. Such, 
basis sets are not 5-function orthonormalizable because the overlap integral between 
the (shielded) spherical Bessel and Neumann functions are not 5-function representable. 
Consequently, the U-matrix in a basis containing such functions is non-Hermitian. This 
is a well recognized result from Kohn calculations. As noted earlier, the A-matrix in 
this basis is also non-Hermitian. 

Through the use of Green’s function, the solution of the Lippman-Schwinger equa- 
tion automatically has the correct boundary condition. Thus the basis for the trial 
functions in a Schwinger calculation need not contain both spherical Bessel and Neu- 
mann functions. It can be a pure L 2 basis or a L 2 basis plus 5-function orthonormal- 
izable continuum functions. In this representation, the fT-matrix is Hermitian. Even 
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if the initial trial basis is not properly orthonormalizable, Appendix A shows that the 
wave function from the next iteration in the Lippmann-Schwinger equation is again 
6-function orthonormalizable. Thus the iterative procedure can be used to ensure the 
Hermiticity property. 


It remains to be shown that the operator A is Hermitian and the A-matrix, rep- 
resented by a properly orthonormalizable basis, is also Hermitian. We recognize that 
H, V, VGpV , and P are all Hermitian operators. On the other hand, the product of 
two Hermitian operators need not be Hermitian. For example, PH is not Hermitian 
because the operator P, defined in the iV-electron space, and H , defined in the N + 1- 
electron space, do not commute. 32 Thus the Hermiticity of A depends on the proof that 
both PH + HP and PV + VP are Hermitian. Consider two arbitrary N + 1-electron 
functions / and g, which can be properly orthonormalized. 


(f\PH\g) = Y, / Is)*. (15) 

rn 

The subscript N under the bra-ket sign denotes integration over the IV-electron space. 

( 9 \mf)' = J df N+i ^(g\$ m ) N (<t> m \H\f) N ^ f = (f\HP\g). 


It is readily seen that 


(g\PH + HP\f)' = ( f\PH + HP\g). (16) 

Thus P H + HP is Hermitian. The proof for PV -1- VP is similar. We conclude that 
A is a Hermitian operator independent of the choice of a as long as it is real. The A- 
matrix is also Hermitian in a basis which is ^-function orthonormalizable. Under such 
circumstances, Eq. (12) is not necessary to maintain the Hermiticity of the A-matrix. 

The present result demonstrates that, either by using a basis containing L 2 func- 
tions plus ^-function orthonormalizable continuum functions or purely L 2 functions, 
or by using a iterative procedure to insure the orthonormalizability of the wave func- 
tion, the Hermiticity of A is automatically satisfied. Then the Hermiticity require- 
ment is insufficient to determine the parameter a, and a becomes a free parameter 
which can be chosen variationally so as to optimize the scattering calculation. This 
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is the approach used in the present study. Alternatively, one may include both the 
regular and irregular free-particle functions in the basis. Without iterating the pro- 
jected Lippmann-Schwinger equation, this results in a non-Hermitian A-matrix. Either 
Eq. (12) or Eq. (14) can be used guarantee its Hermiticity. Finally, as in the complex 
Kohn approach, one may allow the A-matrix to be non-Hermitian. This will then allow 
the parameter a to be varied. 

The use of a pure L 2 basis requires further consideration. In the conventional form 
of the Schwinger variational principle, the Schwinger functional is given by 

,mn = + (*tr ) i' , is.) - - v'c l+ V|4<,+>)J. 

All matrix elements in the above expression involve the operator V. If V is short range, 
then the long tail in the continuum portion of the trial function does not contribute to 
jmn and it is justified to use a pure L 2 basis in this type of variational principle. The 
Schwinger functional in the SMC method, Eq. (11), includes matrix elements of the 
kinetic energy operator. Takatsuka and McKoy 11 pointed out that, by setting N + 1, 
the open channel matrix elements of the kinetic energy operator do not contribute to 
fmn (see Eq. (13)), and a pure X 2 basis again can be employed. In the present approach, 
a is treated as a variational quantity and Eq. (13) no longer holds. However, Appendix 
B shows that, for problems involving short range potentials, the matrix element in 
fmn where the kinetic energy operator operates on purely continuum functions exactly 
cancels the difference between the total energy and the target energy.. The remaining 
matrix elements of the kinetic energy operator always include a bound function in 
the integrand and the continuum character of the trial basis is again unimportant. 
Consequently, even if a ^ N + 1 it is justified to use a pure L 2 basis for this type of 
variational principle. 
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III. //-electron and iV+l-electron projection operators 


Before deriving a variational equation for a, it is worthwhile to consider first the 
difference bewteen an //-electron and N + 1-electron projection operator. The in- 
troduction of a originates from this difference. The Feshbach projection operator, 28 
which operates on the total .Y+l-electron system, is defined without any extraneous 
parameters. 

p N+i = ^ |$ m (l>2, . . . , jY, JV + i)){^ m (l, 2, . . . , iV\ JV + 1)1 (17) 

m 

Here m sums over the open channel functions. If is an open channel function, 


P N +^i = $i. ( 18 ) 

On the other hand, when the //-electron operator P } defined in Eq. (2), operates on an 
//-fl-electron wave function, the antisymmetric nature of fermion wave functions causes 
the result to be more complicated than Eq. (18). Consider the simplest case when the 
target wave function is expressible in terms of a closed shell, single determinantal 
wave function, 

Hh 2, . • . , N) = A | l)*,/?(2) . . . <t> N/2 a(N - 1 )<f>N/20(N )} . (19) 

Here a and f3 denote spin functions. If the wave function of the continuum electron, 
fi(N + 1), is orthogonal to all the target orbitals, then the S z = 1/2 component of 'fj 
becomes 

*i(l, 2, • • • , N,N + 1) = A |^a(l)^/3(2) . . . 4> N/2 l3(N)f i a(N + 1) j , (20) 

We have 

Ptf i (l,2,...,W,W + l) = - 7 =l==$ i (l,2,...,W)/ i a(iV + l), (21) 

V IV + 1 

and 

{9j\HP\V,) = (22) 

A comparison, of Eqs. (18) and (21) clearly shows the origin of the projection parameter. 
When an W-electron projection operator is applied to an W+l-electron wave function, 
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a projection parameter needs to be introduced to account for the difference between N- 
electron and iV-f 1-electron antisymmetrization. Using Eq. (22) it is straightforward to 
derive the result of Takatsuka and McKoy 11 for elastic scattering in the static-exchange 
approximation using a closed shell SCF target function or an open shell target function 
with a simple spin coupling scheme. 

Eq. (21) is not applicable when the orthogonality constraint between the continu um 
electron wave function and the target wave function is relaxed, and/or when the target 
is described by more sophisticated wave functions. Consider the simple case of X 1 !!* 
— b 3 E+ excitation of H 2 by electron impact. Let the wave function for the state 

be represented by 

$Z = A{l<r a a(l)l<r g /3{2)} (23a) 

and the three spin components of the b 3 £+ state by 

= *4{lcr g Q:(l)l<r u a(2)} 

*w> = 4{lff,(l)lff.(2)-^(a(l)fl(2) +/J(l)a(2))} 

* 1-1 = 

It is well established that in calculating this transition, correlation terms of the 


form 




T a = .A{la s a(l)la u a(2)la u /?(3)}, 

(24a) 

and 




'f'b = ^4{lcr ff a(l)lcr s /3(2)lcr u a(3)}, 

(246) 


should be included. These are usually called ‘penetration’ or ‘recorrelation terms’ in 
other methods because they relax the enforced orthogonality between the continuum 
and target functions. Let the projection operator P be generated from the four target 
functions in Eqs. (23a-d). We then have 

P*a = -J=$ 61 (l,2)l<r u /?(3) - -^$ M (l,2)l<r u a(3). (25) 

This simple example illustrates the fact that, when the SMC method is applied to prob- 
lems other than static-exchange scattering with simple target functions, the rationale 
of setting a = N+ 1, based on Eq. (21), is not valid. 


(236) 

(23c) 

(23d) 
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While the use of an iV-electron projection operator leads to the introduction of a 
projection parameter, it does have significant advantages in retaining a simple repre- 
sentation of the projected Green’s function, Eq. (6), Note that the use of the Feshbach 
operator results in a much more complicated expression for the projected Green’s func- 
tion because of the nonorthogonality between the continuum and target functions. 
Furthermore, an W-electron projection operator obviates the need for orthogonality 
constraints between the continuum and target functions. The method has built-in 
nonorthogonality because the plane waves in the homogeneous solution 5 n and Green’s 
function G p 1 ^ are not orthogonalized to the target functions. Thus T a and in 
Eqs. (24a, b) act as correlation terms in the SMC method whereas in other methods 
they play dual roles of correlation and penetration functions. In addition, the com- 
plementary space in the SMC method is directly constructed from 1-aP. In contrast, 
some applications of the optical potential method construct the complementary space 
separately from the P- space. These advantages of the SMC method cannot be over 
emphasized. 
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IV. Variational determination of a 


Sec. II shows that, by choosing a properly orthonormalizable basis, or by iterating 
the projected Lippmann-Schwinger equation, a becomes a free parameter. An optimal 
choice of a should then be derived variationally. To demonstrate that a indeed acts as 
a free parameter, consider using as a trial function the converged static-exchange wave 
function obtained using another method. 

*SE = A$ 0 (l,...,N)f SE (N + l), 

{hsE — -k 2 )fsE — 0, 

and 

se) = 0 . 

Since 'i se contains only open channel functions, 

(S , Ss|J , J3 r |*SB) - (*sb\HP\* S e) = 0. (27c) 

We have 

(T 5 £:|A|^s£;) = (tfssl^PV + VP) - FGpVI'f'si?}. (28) 

The resulting A-matrix is a-independent. Also, Eqs. (27b, c) are equivalent to Eq. (14) 
and the Hermitian representation of A is guaranteed. 

For the second example, consider the trial function obtained from a close coupling 
calculation (which may include closed channel functions), 

Cc\H\$ cc) - 0. (29) 

(^cclAl^cc) = ('Fcc|^(PV + VP) - VGpV - i( PH + IIP)\^ C c). (30) 

The a-dependence disappears. Since Eq. (14) may not be satisfied, whether A in this 
representation is Hermitian or not depends on the original basis used for T cc . 

As a third examples, consider calculation of the lowest pole of H“. In this case 
the open channel space is null. 

P = 0. 


(26) 

(27a) 

(276) 
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The SMC equation in Eq. (9) is reduced to 

-H*n = 0. 
a 

Here ~ is just overall multiplicative factor which can be discarded. We then retrieve 
the Schrodinger equation, independent of a. 

The above examples demonstrate that the choice a = N + 1 is not unique. Except 
in the bound state case where a drops out completely, a variational determination of a 
is generally preferable. Here we consider the simultaneous optimizations of a together 
with the trial scattering functions. Let and ^ be expanded in a set of trial 

functions f r) 



I®1 +) > = £m£.)|/,), 

r 

(31) 

and 

3 

(32) 

The Schwinger functional is now written as 


/ = 

/mn 2 tt 

r 3 r s 

. (33) 

The requirement that / mn be stationary with respect to first order variations in 6 r , c 5 , 
and a, 5/ mn /3t r = 0, 9/ mn /3c 5 = 0, df mn /da = 0, gives the following relationships, 


(S m \V\f r ) -^2c,(f,\A W \fr) =0, 

3 

(34) 


r 

(35) 

and 


(36) 


r s 


The solutions for 6 r , c s are 33 

6 r = £D r ,(/,|F|S„), (37a) 
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(374) 


c, = £(Sm|V|/r)X>.. 

r 

with 

= </^|A‘+)|/,). (38) 

Notice that in the expression for (£? -1 ) rj , the parameter a appearing in is still 

undetermined. Substitute Eqs. (37a, b) into Eq. (36), we get 

EEEE( 2 3 4 5 »^i« B *-</-w«A.,.(/,.ms.>> = o. (39) 

r s t t* 

Eq. (39) is an algebraic equation for a. Using the value of a determined from Eq. (39) 
to evaluate D, we obtain a variationally stable expression for / m „, 

/ mn = -^J2(S m \V\f r )D r .(f,\V\S n ). (40) 

r s 

This expression of / mn is identical with those obained by Takatsuka and McKoy, 11,12 
except for the value of a used. 

Equation (39) provides variationally optimal values of a. It is a high-order algebraic 
equation, the order being proportional to the number of trial functions used. A number 
of useful properties of the solution can be deduced by studying the structure of Eq. (39): 

(1) If the trial functions are the variational solutions of the Schrodinger equation, 

(fs\H\fr) = 0, 

then Eq. (39) will be satisfied independent of the choice of a. This is a physically 
intuitive result and has already been illustrated by the first two examples cited in the 
beginning of this section. 

(2) The solution depends on the initial and final channels under study. Thus in a 
multichannel calculation, different a’s will be used for elastic, superelastic, and inelastic 
scattering. 

(3) For a given transition, the value of a depends on the trial functions used. 

(4) If a complete basis is used, the result should be independent of the choice of a. 

This comes from property (1) and the invariance of Eq. (39) to a unitary transformation. 
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It is important to note that the term basis denotes both the basis set used to expand 
the one-electron functions, e.g., a Gaussian basis set, and the jV+l-electron basis that 
describes the multichannel nature of the wave function, i.e., the open and closed channel 
functions. This property is found to be very useful in testing ihe completeness of the 
basis used. 

(5) The solution of Eq. (39) depends on both the energy and direction of the inci- 
dent and scattered electrons. However, if we consider the following operator equation, 

E E E E v \h)D„U.\H\U)D rt .(U\V = 0 , (41) 

r j t V 

the a s which satisfy Eq. (41) will automatically satisfy Eq. (39). They represent 
a subset of the solutions of Eq. (39) which are independent of the direction of the 
incident and scattered electrons. Also, the a’s determined in this manner satisfy both 
dfmn/da = 0 and dcr^n/da = 0, with a mn denoting the integral cross section. 

In general, the dimensionality of Eqs. (39) and (41) makes a direct solution difficult 
in practical problems. Alternatively, a search procedure can be used. Starting with 
a trial value of a, the matrix elements of can be determined and the scattering 

amplitude / mn calculated using Eq. (40). The parameter a is then varied to search 
for optimal values. Also, we combine the search for the stability of f mn with a test 
of the adequacy of the basis (property 4) by looking for solutions where / mn or cr mn 
is invariant to a. Numerical examples using this technique are presented in the next 
section. 
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V. Numerical Examples 


As an illustration, we studied e-CO elastic scattering in the static-exchange ap- 
proximation. The calculation uses an uncontracted {Cl0s6pld\O10s6pld\C M2s3pld} 
basis with CM denoting the center of mass. The C and 0 basis start with the 9.s5p 
basis of Huzinaga 4 supplemented with one diffuse s and one diffuse p functions, with 
exponents 0.04736 and 0.03654 for C and 0.08619 and 0.06368 for 0. The d exponents 
for C and 0 are the recommended values of Dunning. 35 The exponents for the center 
of mass functions are selected from the table of Nestman and Peyerimhof. 36 They are 
0.032426 and 0.009211 for s, 0.041975, 0.01997, and 0.009467 for p, and 0.04349 for d. 
All calculations were carried out at the internuclear distance of 2.132 a 0 . 

Sec. IV shows that the SMC static-exchange result should be independent of a if 
the one-electron basis has converged. Fig. la and b present the integral cross sections 
at 3.75 eV incident energy as a function of a, for a = 10 - 20 in steps of 0.1. (The value 
of a is 15 if it is set to N + 1). Of 100 points calculated, 81 fall within 5% of each other, 
as depicted in Fig. la. Fig. lb shows the full calculation where a number of significant 
oscillations in the cross sections are shown. These instabilities are a direct consequence 
of the incompleteness of the basis used. The result of the a variation means that the 
basis set is probably adequate if 95% accuracy in the cross section is acceptable. 

The source of the sharp peaks and dips in Fig. lb is analogous to the pseudo- 
resonances in the Schwinger method reported by Apagyi et al. 3T and analyzed by 
Weatherford et al. 38 and by Winstead and McKoy. 39 When the basis for expanding 
the wave function is complete, the matrices of i T, PH and HP all have zero eigenval- 
ues and the eigenvalues of in this representation do not go to zero unless a real 
resonance exists. On the other hand, the matrix of H represented by an incomplete 
basis generally does not have zero eigenvalues. Thus a particular choice of a may cause 
a near cancelation between H/a and the remaining terms in A ^ such that A has an 
eigenvalue close to zero. Similarly, a dip indicates a corresponding maximum eigenvalue 
^ Notice the sharpness of these structures, a characteristic of pseudo-resonances. 
Since dcr mn / da 0 is also satisfied at these structures, we consider a choice of a based 
on the flatness of the cross section vs a curve to be superior to the simple condition 
do-mn/da = 0 . 
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As a second example, we studied the effect of the variation of a on the e-H^, 
— ► b 3 E+ excitation cross sections using a two-state approximation. This is a well 
studied transition; calculations using the SMC, 15,20 R-matrix, 40,41 linear algebraic, 42 ’ 43 
and S-matrix Kohn 6 ’ 44 methods have been reported. These calculations showed that 
the correlation/penetration terms such as those given in Eqs. (24a, b) are important in 
describing the scattering process. 

Two Gaussian basis sets were used. The 64-basis, consisting of 64 uncontracted 
Gaussian functions, is a {H6s6p|CM4s4p} set that was used in the original calculations 
of Lima et al 15,20 and the cross sections are in good agreement with experiment. 45 ” 47 
The 70-basis, composed of 70 uncontracted Gaussians, is a {H8s7p|CM3s3p} set which 
employs the s and p functions from van Duijneveldt’s 48 8s6p4d hydrogen basis, sup- 
plemented with one p function (exponent 0.03125) at the nuclear centers, and three 
s functions (exponents 0.082, 0.021, 0.0052) and three p functions (exponents 0.114, 
0.0226, 0.0045) at the center of mass. The hydrogen basis of van Duijneveldt, in the 
form of atomic natural orbitals, has been used frequently in high quality quantum 
chemical calculations of hydrides and their ions. 49 Based on previous experience, we 
expect the 70-basis to be superior to the 64-basis. The present study employs the same 
type of target functions as previous SMC calculations on this transition, 15,20 i.e., an 
SCF function for the ground state and an IVO function 50 for the b state, so that the 
a-dependence of the cross section would not be attributed to the difference in the target 
wave functions. All calculations were carried out at a bond length of 1.4 a 0 using the 
AMES SMC code. We determined the integral cross sections as a function of both a 
and the electron energy. The range of a used is 2.0 - 4.0 in steps of 0.05. Also note that 
the results of the 64-basis calculations are slightly different from Lima et al. 20 because 
a better representation of the V G^V term is used. 

Analogous to the CO static-exchange result, the e - elastic cross sections versus 
a curve also show regions of stability, with sharp structures between stable regions. Figs. 
2a and 2b depict the elastic cross sections at 13 eV, calculated using the 70-basis and 
64-basis, respectively. Only the data points within 10% of each other are depicted. For 
the 70-basis, 72.5% of the data points satisfy this criterion, versus 45% for the 64-basis. 
From this point of view, the 70-basis is superior, even though the cross section averaged 
over the converged data points, shown as a dashed line in the figures, is in agreement 
with each other, 23.4 au for the 70-basis and 22.8 au for the 64-basis. Even for the 
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70-basis, the convergence is inferior to the CO result. We consider this a consequence of 
the inadequacy of the N + 1-electron basis used in the two-channel calculation, and not 
due to an inferior one-electron basis. Note that the N + 1 -electron basis (multichannel) 
does not play a role in the static-exchange approximation used in the CO calculation. 

Figures 3a and 3b depict the inelastic cross sections from the 70-basis and 64- 
basis calculations where the data points are within 10% of each other. Among the 
calculated values, 45% and 50%, respectively, fall in this range. Thus for the 70-basis 
the convergence of the inelastic cross sections is significantly inferior to the elastic case 
whereas the quality of the 64-basis cross sections remain roughly the same. The 70-basis 
result serves as an indication that an optimal basis for elastic and inelastic calculations 
may be different. 

Note that for the 64-basis, the regions of stability for both elastic and inelastic cross 
sections include the point at a=3.0, corresponding to Takatsuka and McKoy’s choice of 
^ ~b 1 1 whereas the cross section at this a is not in the set of converged points for the 
70-basis. Fig. 4 compares the inelastic cross section as a function of electron energy, for 
the 64-basis and 70-basis at a = 3.0, and the 70-basis cross sections determined using 
the average of the converged cross sections at each energy. Experimental data 45-47 for 
this transition are also included for comparison. Note that the 64-basis cross sections 
at a 3.0 and the average cross section from the 70-basis calculations are in agree men t 
with experiment and with each other to within the estimated errors, but the 70-basis 
cross section curve at a=3.0 differs significantly in both shape and magnitude from 
the other two curves and from experiment. In particular, the experimental data show 
no evidence of a dip around 16 eV. If a = N -f- 1 is considered a unique choice, we 
have the paradoxical result that a larger and, generally considered better, basis gives 
cross sections that disagree not only with the previous SMC results, but also with 
calculations using other methods and with experiment. However, if we accept a as 
a variational parameter and the optimal cross section using a-stability criteria, then 
the 70-basis results agree with the generally accepted values. This numerical example 
demonstrates the importance of a- variation. 

All previous reported SMC cross sections 13 ” 26 were calculated using a = N+l. 
From our point of view, this is just one of many possible values of a. In general, a 
calculation of the cross section vs a curve is recommended since it provides both a con- 
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vergence test and a ready check for pseudo-resonances. Furthermore, the requirement 
that the cross section curve be flat with respect to a provides a variationally more 
meaningful cross section than calculations using one single value of a. 
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VI. Conclusions 


This study considers the role of the projection parameter in the SMC method and 
develops a variational equation for it. While the multiple solutions in the variational 
equation may at first appear to be a disadvantage, it actually provides a sensitive test of 
the convergence when we couple the variational solution of a with the requirement that, 
for a fully converged calculation, the cross section is invariant with a. The convergence 
requirement applies both to the one-electron and N+l -electron basis. Such a test 
is apparently unique to the SMC method. It is also shown that calculations using an 
incomplete basis, as is the case in all practical calculations, may find pseudo-resonances 
when the choice of a causes the SMC operator to have a near zero eigenvalue. Thus it 
is essential to test the stability of the cross sections with the a-variation. Studies on the 
invariance of the cross section with a using more sophisticated target wave functions 
and more extensive multichannel treatments are under way. 
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Appendix A. Orthonormalizabilty of the iterative solutions of 
the projected Lippmann-Schwinger equation 

Consider the projected Lippmann-Schwinger equation 

p^r(+) = ^(+) +G .(+)y^,(+) ) (4-1) 

where + 1 j s the projected homogeneous solution and satisfies the 6-function orthonor- 
malization condition, 

<*4 +, |v4 +) > = - £,). (A- 2 ) 

In Eq. (A-l), is the wave function represented in the initial trial basis and is 
the wave function after one iteration of the projected Lippmann-Schwinger equation. 
The initial trial basis may be composed of purely L 2 functions or L 2 plus continuum 
functions. If continuum functions are included in the basis, 6-function orthnormaliz- 
ability is not required. Consider the overlap integral 

(P4 (+)| P4 (+ ) ) = ( v ,(+)|p<+>) + (»,‘ + '|G ( p +) (J5 J )V't<t > ) 
+{Gp +, (£i)V''t ( ,t ) lv=i +> ) + (G , r HE l )V^ > \G < ,t > (Ei)v¥+ l ). (A - 3) 

Direct integration gives 

i — (A -4) 

(A- 5) 

and 

( o' + 1 ( £, j F i ’ I G p 1 ( £, j V >1 < +>) 

= tim,-o ^ ^ + — [(Gg~ ) (gi)V*i< ( ' l ' ) |V|*at ) ) - <*ifi|V1G ( p + \E 2 )V<M< + >)]. 

(A-6) 

Combining Eqs. (A-2), (A-4), (A-5), and (A-6) and taking the limit 6 =0, we find 

(p 4 t+)|p$(+) )={ (J 1 _£ J ) + _^i ( ^(+> |K |^<+ ) )_( 4 («| F |p 4 (+ ) ) ]. (jt -T) 

The presence of the interaction potential V, which vanishes at least like r -1 as r — + oo, 
guarantees that the last two matrix elements in Eq. (A-7) do not diverge for neutral, 


20 



bound target functions even if and are represented in a basis including 

continuum functions which are not 5-function normalizable. Eq. (A-7) means that 
\ and hence are 5-function orthonormalizable. Thus the iterative solution 

of the projected Lippmann-Schwinger equation is 5-function orthonormalizable even if 
the initial trial function is not. 
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x*xAknrt)Y Xn „(r 3 ) (5-4) 

The first term on the right hand side of Eq. (B-4) vanishes due to Eq. (B-2). The 
second term involves the short range potential V which does not weigh the long range 
part of s\ m and The last three terms on the rhs of Eq. (B-4) comes from the non- 
orthogonality between the target functions and the spherical Bessel functions. The 
transition one-, two- and three-particle density of target functions <f? m and 4> n are 
denoted by Dmh, Plnn, and Dmi, Pi, P 2 and Pj represent the number of permutation 
operations, and U(f) the nuclear attraction potential. Note that ail target transition 
density functions are short-range. Thus the last three integrals in Eq. (B-4) involving 
short-range integrands in all coordinates. The long tail of the continuum functions s\ m 
a nd -sa„ again does not contribute to these matrix elements. It is straightforwad to 
show that similar conclusions apply to matrix elements of PH and HP. 

In view of the above result, the matrix elements in the Schwinger functional in 
Eq. (11) is insensitive to the continuum nature of the trial function even when a/ff-f 1 
and Eq. (13) does not apply. This justifies the use of pure LP basis in the variational 
principle developed in this paper. 
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Figure Captions 

Fig. la e-CO elastic scattering cross section in the static-exchange approximation versus 
a at 3.75 eV electron energy. This figure depicts all data points where the cross 
sections are within 5% of each other. 

Fig. lb Same as Fig. la except all data points are included. 

Fig. 2a e-H 2 elastic scattering cross sections at 13 eV as a function of a. This figure depicts 
all data points from the 70-basis calculations where the cross sections are within 
10% of each other. The dashed line denotes the average cross section over the data 
points. 

Fig. 2b Same as Fig. 2a for cross sections determined using the 64-basis. 

Fig. 3a e-H 2 — ► b 3 S+ excitation cross sections at 13 eV as a function of a. This 

figure depicts all data points from the 70-basis calculations where the cross sections 
are within 10% of each other. The dashed line denotes the average cross section 
over the data points. 

Fig. 3b Same as Fig. 3a for cross sections determined using the 64-basis. 

Fig. 4 e-H 2 X*E+ — + b 3 E+ excitation cross sections as a function of electron energy. The 
experimental data are from Hall and Andric (Ref. 45), o; Nishimura and Danjo 
(Ref. 46), o; and Khakoo et al. (Ref. 47), A. 
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ABSTRACT 

The Schrodinger equation for the scattering of an electron by a hydrogen molecule 
is solved by converting it to a block-tri-diagonal matrix equation. The discrete mesh of 
points in the angular variable must be made more dense as the potential energy deviates 
from spherical symmetry to maintain a desirable level of accuracy. The present paper 
investigates the nature of this dependence, and the relative effect of third-order and fifth- 
order differentiation rules. It is found that the existence of one large element in the 
K-matnx sets a lower bound on the accuracy of all elements. This bound, and therefore 
the overall accuracy of the numerical approximation, can be substantially improved by the 
use of the higher-order rules. 
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1. Introduction 

Traditional electron scattering theory concerns the results of collisions between an electron 
and a target under the assumption that their interaction depends only on their separation, and is 
independent of their angular orientation. When the target is an atom, this assumption is justified, 
but when the target atoms are bound into molecules, the great simplifications afforded by the 
conservation of angular momentum are no longer available. A numerical analysis must therefore 
be introduced at an earlier stage in the solution of the problem. 

The K-matrix is a concise and computationally simple means of presenting the results of a 

scattering calculation done by the method of partial differential equations. This work will 
investigate the relationship between the mesh size required in a finite-difference approximation to 
the K-matrix and the deviation of the scattering potential from spherical symmetry for various 
levels of accuracy. "Accuracy" is defined as the difference between the numerical approximation 
and a hypothetical continuum-limit K-matrix, which must be inferred from the behavior of the 
approximation. 

The angular accuracy is roughly independent of the radial properties of the potential, so the 
numerical requirements of the polar-angle dependence can be investigated without a detailed model 
of the physical potential. The angular momentum quantum number / provides a convenient basis 
for calculating the K-matrix, but one which is not a basis of eigenvectors. In an eigenvector basis, 
the K-matrix would be diagonal, and the numerical convergence of the various elements would be 
unrelated, with the states of higher rotational energy requiring a finer mesh. Here, by contrast, the 
enror in the numerical approximation is spread over all coupled states. In particular, if one element 
dominates the others, a small relative error in that element will imply a large relative error in the 
other elements. The convergence of the largest element is found to be excellent, with accuracy 
about 10-3 using 29 points in the angular mesh. The convergence of the smaller elements is limited 
by this effectively constant bound. Convergence beyond this limit requires vastly larger meshes, 



3 


though a great improvement can be gained by the use of high-order differentiation rules. It is 
found that adding rules of fifth order in the mesh spacing requires little additional programming, 
and virtually no extra computer time. 

2. The Method 

The homonuclear diatomic molecule is the simplest system for which non-spherical terms 
contribute to the scattering potential, so it will serve as the physical anchor for an otherwise 
mathematical study. When the energy of the incident electron is between 1 and 10 electron volts, 
the non-spherical structure of the molecule is particularly important. In such a situation we use the 
axial symmetry of the molecule to reduce the Schrodinger equation to a two-dimensional partial 
differential equation. We assume for simplicity that the collision does not excite vibrational states 
of the molecule. Rotational excitations may be dealt with at the cross-section stage of the 
calculation in the usual adiabatic-nuclei manner [1], 

The target molecule may be assumed to be oriented in a particular direction in space, so we 
use the axis of symmetry of the molecule as tire z-axis of our coordinate system (Figure 1), and 

define a modified wave function u which is related to the true wave function d> for the projectile by 

0(r, d, 4>) = r-Ju(r, d) eimt. 

The equation for u in atomic units is 

d 2 u f du If m 2 ) 

+ rMcotd + \ + { p 2 _ 2 V(r,d) U =0. (2) 

3r2 l 3d J 1 ' r 2 sind J 

where V(r, d) is the scattering potential, and p is the momentum of the projectile. A finite 
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difference approximation to this equation is solved on the domain 0 < r < R max , 0 £ i3 £ rc/2. The 
outer limit R ma x is chosen to be sufficiently large that the scattered particle is effectively free at that 
distance. The range of $ has been cut in half to take advantage of the symmetry of the molecule 
under reflections in the x-y plane. 

The boundary condition at r = 0 requires that u= 0. At t!M), thbre are two possibilities: if m, the 
z-cpmponent of the orbital angular momentum, is nonzero, then n=0; otherwise 8u/Bd=0. At 

i3=7t/2, the parity of the system dictates either that u = 0 or 3u/9d=0. At r=R max , we impose the 
requirement that the incoming electron be in a state of well-defined orbital angular momentum, so 

u(Rmax. t3)=Pl m (C0S l3). (3) 

The parity of the associated Legendre function, (-l)'+ m , determines the choice of boundary 

condition at i3=k/ 2. Though this parity is conserved by the scattered wave function, the value of / 

for the incident wave is not. Outgoing waves therefore contain a superposition of /, /±2, /±4, and 
so forth. 

Once the wave function has been found, information about the scattering is extracted by 
matching it to the known asymptotic form 

= X [Au'ji'fpr) + B]ryi(pr)] Pj m (cos i3), (4) 

r 

where / and m are the index and order of the Legendre function used in the boundary condition at 
large r, andy and y are spherical Bessel functions. The coefficients A and B form matrices on the 
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space of allowed values of orbital angular momentum. 

The orthogonality of the Legendre functions may be used to project out a particular value of 
When this is done for two different large values of r, the linear independence of the Bessel 
functions allows extraction of the unknown coefficients A and B. This process must be done for 
each allowed value of /'; then, the whole solution process repeated, allowing / to take on each 
permitted value. The result is two matrices describing the transformation of an angular momentum 
eigenstate into a scattered wave. In principle, the dimension of the matrices is infinite, but in 

practice the finite range of the non-spherical potentials provides a cutoff. For an H 2 molecule, 

with a projectile energy below 10 eV, values of / less than 5 are all that are required. Hence, the 
matrices are at most 3x3. The K-matrix is defined by the equation 

^lm = X [}\ipr) + Kj]'yi'(pr)] P r m(cos i3), ( 5 ) 

1 ’ 

so K = A-1B, and the influence of the artificial imposition of an initial orbital state is removed. The 
diagonal elements of K are the tangents of the elastic phase shifts, and the off-diagonal elements 
describe the coupling of states of different orbital angular momentum by the potential. 

3. The Potentials 

The potential energy function V in equation (2) represents the electromagnetic interaction of 
the projectile with the target molecule, and is the source of the breaking of spherical symmetry. In 
principle, this interaction is highly complicated. A molecule is by no means a static, rigid body in 
the presence of an incoming electron, so a potential model is not strictly correct, but the practical 
usefulness of potentials in calculations is great enough to excuse the gross simplification of the 
scattering problem they imply. The potentials involved in scattering from hydrogen molecules are 
approximations to three effects: the static electric moments of the molecule; the polarization of the 
molecule by the electric field of the projectile; and the substitution of an electron from the molecule 
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for the projectile in the outgoing state, known as "exchange". Toy models have been adopted 
which mimic to some extent the radial structure of the physical potentials, and have angular 
properties which are expandable in the low-order Legendre polynomials of even parity. 

Because the hydrogen molecule is composed of identical atoms, the lowest-order multipole 
field seen by the projectile is the quadrupole term 

. V Q (r, i3)= Q r-3 P 2 (cos J), rg) 


where Q-.49 in atomic units. This form is valid for values of r which are outside the molecule for 
all A At shorter ranges, the quadrupole potential is cut off by the usual exponential term 

(l-exp[-(r/a)6J) as a multiplier. The adjustable parameter a is set to 1, since no effort is made to fit 
experimental data. The static potential inside the molecule is represented by a Yukawa function of 
unit range and arbitrary strength, which is adjusted as necessary to give the correct order of 
magnitude for the K-matrix. 

Polarization of the target is well known to include terms of both spherical and P 2 character 
[1]. The form used here is 

Vp(r, i3) = H(l-exp[-r6]) [a 0 + a 2 P 2 (cos fl)]/2, ( 7 ) 

with a o=5.5, and a 2 =\A. 

The exchange potential is a crude form of the Free Electron Gas model. The basic idea here 
is that the electrons in the target may be approximated by a degenerate Fermi gas whose density 

vanes with position [1]. The exchange potential is proportional to the Fermi momentum at any 
point in space, giving 


V x (r) = -(3/27t) [3rc 2 B(r)]i/3. 


( 8 ) 
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The electron density B(r) is approximated by a simple sum, 

B(r, i3) = (4n)-ie-r[l + 1.3 P 2 (cos i3)], (9) 

which has contours of constant density which resemble those of the sum of two hydrogen atoms 
separated by .7 angstroms. The resulting potential can be approximated by a linear combination of 
Legendre polynomials, which is 

V x (r,fl) = -.8 e-r/3 [.97 + .43 P 2 (cos fl) + .1 P 4 (cos fl)J. (10) 

The "realistic" potential is a sum of these three terms. 

The long-range quadrupole potential is the most important term in real scattering problems, 
but it does not vanish until a distance of 14 Bohr radii from the molecule's center is reached. This 
was determined a posteriori, by searching for a region in which the K-matrix does not depend on 
small changes in R ma x-) Such a long range requires a large number of radial points, and 
consequently a long running time for the computer code. Although the convergence properties of 
the angular approximation are roughly independent of the quality of the radial approximation, this 
independence is only valid when the radial mesh is on the order of magnitude of the mesh which 
gives decent radial convergence. Therefore, the quadrupole term is generally omitted from this 
work, in the interest of a more efficient use of computer time. The P 2 terms in the polarization and 

exchange potentials are in any case expected to yield the same requirements on the angular mesh. 


4. Numerical Analysis 

The usual discrete approximation to equation (2) is obtained by replacing the radial variable with 
a set of N r points, and the angle with a set of N 0 points. The wave function therefore becomes a 
set of (real) values which solve the equation [2] 
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[ u i+l,j ~ 2Ujj + Uj.j j ] + 

r i' 2 [ h o _1 COt (tljj +1 - Ujj) + h 0 -2( Ujj +1 - 2Ujj + Uy.i) ] + 

[p2-2Vij- m2 / ri 2sin20j] ujj = 0. (11) 

Here, My-ufa, t)j), with 0 < j ^ N r and 0 <j <L N 0 . The spacings between the radial and angular 

points are h r and h Q , respectively. (The values of h t and h 0 may in principle depend on i and j.) 
Equation (1 1), together with the boundary conditions, gives a system of (N r +2)(N 0 +2) coupled 
linear equations for the «y, which may be solved by conventional matrix methods [3]. 

The choice of the points r ; and i3j presents an unavoidable conflict of interest. There are 

ways to improve the accuracy of a numerical solution of a differential equation by judicious choice 
of abscissae, just as there are ways to improve the numerical quadrature by which the elements of 
the K-matrix are projected out of the wave function. Unfortunately, the ways are mutually 
inconsistent. Numerical quadrature is optimized by choosing points so that successive errors 
alternate in sign, which is the worst possible choice of points for differentiation. Therefore, a 
compromise of equally-spaced abscissae will be used, and optimization achieved through 
higher-order integration and differentiation rules. 

When the Schrodinger equation is solved on a two-dimensional lattice, the resulting matrix 
is block tridiagonal. This property, which makes the matrix easier to decompose into lower and 
upper triangular matrices, depends on one of the variables having a simple three-point derivative. 
The other may be as complicated as the problem demands. In this work, the radial accuracy is of 

secondary interest, so the i3-derivative part of equation (11), which is in braces, will be explored 
with rules of greater accuracy (Table 25.2 in reference [4]). 
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For a point in the interior of the domain of approximation of i3, labelled by "j", rules of order 
h 0 5 are 

dm = (l/DO { 2fj_2 — 16fj_i + 16fj+, - 2fj + 2 } , (12) 

a2f/3i32= (1/D 2 ){— fj- 2 + 16fj.i - 30fj + i6fj + i - 2fj +2 } , (13) 

where Di=24/i 0 and D 2 =12/i 0 2 . (To simplify the formulae, u\ j has been replaced by fj.) For 
points next to the boundary, asymmetrical rules must be used: 

3f/3i3 = (1/DO { — 6fj_i - 20fj+ 36fj+i - 12f j+2 + 2f j+3 } ; (14) 

3 2 f/3i3 2 = (1/D 2 ) { 1 lfj.i - 20fj + 6f J+ i +4fj + 2-f j+3 }. (15) 

These rules are used for points next to the boundary at i3=0. Next to i3=7c/ 2, the subscript "j+3" is 

replaced by ‘’j-3”, "j+2" by "j-2" and so forth, and there is an overall minus sign in the first 
derivative. 

At the boundary, if u= 0 is the boundary condition, then it goes directly into the matrix. If 
Neumann boundary conditions apply, the process is more involved. First, terms in the 

Schrodinger equation proportional to 3f/3i3 are dropped. Then, equation (14) is used with j=0 to 
get the relation 

-6^1-20^ + 36^- 12f 2 + 2f 3 = 0 (16) 

which is solved for f_i. This value is then used in equation (15) (with j=0 again), to derive an 
expression for 3 2 u/30 2 at J=0 which is accurate to fifth order in h c : 


02f/ai32 = (1/D 2 ) {— (170/3)f 0 - 72f, - 18f 2 - (8/3)f 3 } . 


(17) 
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When Equation (17) is used in place of the second derivative term in (1 1), the result is an 

approximation which significantly improves the numerical accuracy, compared to the three-point 
rules. As a test, when a spherical potential is assumed, and the K-matrix calculated, the off- 

diagonal matrix elements should vanish. In fact, for a given number of 0 points, they are at most 

one fourth of those calculated with the three- point rules, and of order 10-3 when the largest 
diagonal element is about 1. 

5. Results 

Convergence of the K-matrix is determined in two ways. The first takes advantage of the 
invariance of the K-matrix under time reversal to require that the K-matrix be symmetric. The 

second is that the value of an element of the K-matrix change with changes in the number of fr 

points in a manner which points clearly to an asymptotic limit. 

The matrices are not stricdy symmetrical. When the scattering potential is spherical, the 
wave function should have the angular properties of the large-R boundary condition, all the way to 
the origin. This is true only when the boundary condition is 1=0. When higher waves are used, 
each contains a contaminant of lower / values. This contamination decreases with each step in r 
away from the origin, reaching an eventual limit of about .002 times the amplitude of the correct 
wave. Therefore, any number less than .002 in the K-matrix is indistinguishable from zero. This 

seems to be an artifact of the solution algorithm, reducible only by improving the i3-derivative rule. 

(When a three-point rule is used, the corresponding limit is .008.) 

The second criterion makes an accurate graphical depiction of the convergence difficult. 

The convergence of almost all elements is monotonic, so the choice of a specific limit is arbitrary, 
yet it has a great impact on the perceived quality of the convergence. For this reason, the first 
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derivative of with respect to N 0 is plotted in Figures 2-7. 

To best isolate the effects of various non-spherical potentials, the K-matrices used to 
produce figures 2-7 are derived from several potentials. First, a spherical Yukawa potential was 

used, then the same Yukawa potential multiplied by the Legendre polynomial P 2 (cos 0), then half 

of each. The length scale of the potential was 3, and its strength was set to -0.4. The momentum 
of the incoming particle was .8, corresponding to a kinetic energy of 8.7 eV. All these parameters 
were chosen to give significant K-matrix elements up to 1 = 2 . An exception is the P 2 potential in 
Figure 5. The 1=2 element was very small, so its relative error was extremely large, and not 
indicative of the numerical properties of the algorithm. The fourth curve on the graphs is the 
"realistic" potential described above, with a Yukawa of unit range in place of the long-range 
quadrupole. 

Figures 2 and 3 show the slope versus N 0 of the K-matrix element for transitions between 
spherical waves. Figure 2 is for a three-point rule, and Figure 3 is for the five-point rules. The 
advantage of the five-point rule is most obvious at small numbers of theta points. As the mesh 
becomes finer, the relative advantage of the five-point rule diminishes, but it is still a factor of two 
better. Figures 4 and 5 are the same, for the D-wave to D-wave reaction. 

Figures 6 and 7 are different. They show the difference between the 1=0 to 2 matrix 
element and the 2 to 0 element, and thus provide a test of time-reversal. The difference is 
normalized to the value of K 02 . usually the smaller of the two. In judging the size of the relative 
error, it should be borne in mind that there is always a diagonal element of the K-matrix which 
dominates this term by an order of magnitude. The uncertainty in this dominant element is the 
primary source of this error. In this place, the five-point rule shows its particular virtue. The 
five-point rule shows a more symmetrical K-matrix with 1 1 points in the mesh than the three-point 
rule with 30. 
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6. Conclusion 

The addition of five-point approximations to the angular derivatives requires only minor changes 
in the computer code for solving partial differential equations. It preserves the block-tridiagonal 
character of the equivalent linear system, so it requires no extra running time, and it yields 
significant improvement in the accuracy of the numerical solution. 

The chief drawback of the partial-differential-equation method for non-spherical scattering 
problems is that the coupling between states makes the relative error in the small elements of the 
K-matrix large, since it is a small fraction of a large number. The largest element of the K-matrix 
thus creates a "floor", a lower bound on the possible error in the other elements. 

The contamination of the wave function by lower partial waves than the input boundary 
condition is the largest single source of error. It is an artifact of the «= 0 boundary condition at the 
origin. The natural tendency of this algorithm is to assume the smoothest possible wave function 
near r=0, and it takes several steps away from the boundary before the large-r boundary condition 
can assert its dominance. A more sophisticated treatment of the radial derivatives might well be 
rewarded by a reduction of this defect. 

Though this work has confined itself to the angular behavior of the wave function, the great 
improvements afforded by the higher-order approximations to the derivatives might have 
counterparts in the radial direction. This would destroy the tridiagonality of the matrix, but in 
applications where the size of the computer's memory is less of a constraint than the cost of 
processor time, it may be worthwhile. 
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Figure 2. The relative change in the K-matrix element for Srwave to S-wave transitions, calculated 
using three-point derivative rules. 
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Figure 3. Same as Figure 2, but with five-point rules. 
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Figure 4. The relative change in the K-matrix element forD 
calculated using three-point derivative rules. 
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Figure 5. Same as Figure 4, but with five-point rules. 
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Figure 6. The relative change in the K-raatrix element transitions between S- and D-waves, 
calculated using three-point derivative rules. 



5 



Number of Theta Points 


Figure 7. Same as Figure 6, but with five-point rules. 
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Abstract 

A calculation of e-N 2 scattering in the vicinity of the 2.4 eV (U g ) resonance has been 
completed. The main element of the calculation is a 15 term vibrational close coupling 
expansion, reduced to coupled two-dimensional (2d) partial differential equations (pde’s), 
and solved using the noniterative pde technique. The potential consists of static, exchange, 
and polarization parts; each part has been (previously) derived in a manner appropriate to 
its importance in the scattering equation. Results for the absolute total cross section, both 
in magnitude and shape of the substructure in the resonance region (1.5< k 2 <3.0eV), are 
in excellent accord with experiment. Angular distributions are also calculated and found 
to vary significantly in shape in the immediate vicinity of the center of the resonance 
(2.05<k 2 <2.15eV), indicating the need for differential measurements at a finer grid in 
energy, and therefore, requiring even better energy resolution. Comparison with other 
calculations and discussion of some theoretical aspects are also included. 
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I. INTRODUCTION 


The purpose of this paper is to report a calculation of electron scattering from molec- 
ular nitrogen (e-N 2 ) using the noniterative partial differential equation technique. 1 This 
calculation is confined to the energy region around the 2.4 eV resonance; it consists of a 
15-state vibrational close coupling expansion for the resonant (II 3 ) partial wave combined 
with fixed-(plus adiabatic-)nuclei amplitudes for the other (contributing) partial waves. 
The underlying method is the hybrid theory. 2 We shall concentrate on the angular distri- 
butions for the lower vibrational states (v = 0, 1,2), showing that their shapes change very 
rapidly from one energy to the next as one traverses the energy region in the center of 
the resonance. Thus, despite the fact that several differential scattering experiments have 
been reported (to be discussed below), they do not have the energy resolution, nor have 
they attempted to explore the variation in the angular distributions in this narrow energy 
region in detail. It is one of our aims to motivate such experimental investigations. 

Several important developments in our calculations have been made since the intro- 
duction of the hybrid theory 2 , which we have reported piecemeal since then. Briefly, the 
most salient of them are: the noniterative technique itself 1 , the reduction of the scattering 
equation to 2d form, which was first outlined in Ref. 3, in addition to which— it turns out— - 
the 2d technique could also be applied to the derivation of the polarization potential 4 , as 
well as to a method for exactly including exchange in the static-exchange approximation. 5 
And finally, the diiect (i.e. non-exchange) static potential could be calculated 6 using a 
much better (MCSCF) approximation of the N 2 ground state than the SCF approxima- 
tion, used originally. 2 5 Some interim results using a 10 state vibrational expansion were 
reported in Ref. 7. 


II. THEORETICAL ASPECTS 


We start with an antisymmetrized ansatz for the total wavefunction of the e-N 2 system 

2 , 7 . 
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*<”•> = (2.1) 

i= 1 

Here x, are the coordinates (space and spin) of the i th electron and x ( ') is the collection of 
coordinates of the remaining (fourteen) electrons. The factor (-l) pi is the parity of a cyclic 
pei mutation ( Pi ) of the sequence 1,2,..., 15, thus making completely antisymmetric. 
We first summarize the analysis, mostly of Ref. 5, whereby the basic integro-pde is 
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derived. (Note, however that we are here using Rydberg units as opposed to atomic units, 
used in Refs. [5,6]). Insertion of Eq. 2.1 into the Schrodinger equation, and premulti- 
plication by the target ground state, yields the static exchange approximation. We 

first consider the ground state (E+) to be represented by a single determinant (i.e. SCF 
approximation): 

= det(\a]2o]\al2olZal\*l u Ul u - E+) (2.2) 

Labelling the different orbitals <f> a [a : 1<t 9 , • • • , l7r yu ], we recognize that each cf> a is an 
explicit function of the coordinates of a single electron and an implicit function of the 
intemuclear separation R: 


<t> a = (f> a (r u - ■ ■ -R)( Q (spin) (2.3) 

When one includes the R dependence of (j> a and F^ m \ one derives 7 a 3d pde for F^ n, h 

7 

[-V 2 + H vib (R) + V(r- R) - E sc ] F< m >(r; R) = 2 £ W< m >(r; R)^(r; R), (2.4) 

a=l 

where V(r*;R) is (to begin with) the static potential between the scattered electron and 
the target (N 2 ): 


Ks<a<ic( r i 7?) — (^*^2 |Le — moi|^N 2 ) (2-5) 

and W a are the (static) exchange kernels 

R) = J F<”>(r R) (2.6) 

The energy appearing in Eqs. 2.4 and 2.5 is 

Esc = k 2 v + e„, (2-7) 

and H vib (R ) is the vibrational part of the target Hamiltonian from whose internuclear 
potential, vibrational wave functions (x„) and energies ( e v ) are evaluated. (In this calcula- 
tion, the latter are obtained numerically from the MCSCF potential energy curve reported 
in Ref. 6. 

For the static potential (2.5), we use an MCSCF approximation 6 of An MCSCF 
wave function is a sum of determinants, whereas the right hand side (i.e. the exchange 
part) of the scattering Eq. 2.4 assujTies that is a single determinant (i.e. SCF 
approximation). Thus (2.4) constitutes an admittedly unbalanced approximation, which 


2 



nevertheless seemed well justified 6 - 7 because (a) the exchange terms usually have a quanti- 
tatively smaller effect than the direct terms, and (b) when written as in (2.4), the exchange 
terms can be rigorously reduced to a coupled set of ordinary (i.e. non-integral) pde’s. 6 This 
comes about by expanding F ^ and in vibrational states \v of N 2 


F {m \^R) = Y J Fi™\r) Xv {R) 

v-Q 

N v 

R) = Yl 


(2.8) 


and, using the well known property of the Coulomb potential as the Green’s function of 
the kinetic energy, 


v2 (j>_ ^ /I ) = -4irS(r ' - r), 


(2.9) 


allows the scattering equation (2.4) to be reduced to coupled, but nonintegral pde’s 7 : 


N v 


[V 2 + kl] Fj»>( r) = E 


u y = 0 


V vv .(r)F ( v : n \f) - <t> ( v a J(r)W ( v r' a )(r) 


< c 


a=l 


N „ 


V 2 H 2 <”’’“)(0 = - 4 s ^ *< 2 ( r ) f £">( r ). 


( 2 . 10 ) 


u '=0 


The double indexed quantities in (2.10) are vibrational matrix elements of the unsub- 

scripted quantities. Thus the matrix elements in Eqs. 2.10 are given by the following 
expressions: 


' / '»'(0=( x „(B)|F(f;f?)| x „.(J{)), 

= (x«(-S)IM?; rt)| xv(tf)), 




= E^tOI^I^')) 

v f =0 \ T V I 


( 2 . 11 ) 


Eqs. 2.10 can be further reduced to 2d pde’s by exploiting the cylindrical symmetry of the 
various functions. As derived in Ref. 5, they take the form 


+ kl] = E - i E 

v' = 0 L r 

A ( m - = ~~Y ^15U)A m) U). 


Ot— 1 


( 2 . 12 ) 


u — 0 
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here z — (r, 0) and A (m) is the 2d Laplacian 


A , v d 2 1 


d 2 1 

~ ! 2 + 


<90 


5 9 m 2 

Q ? 77/7 — 

sinO 86 dO sin 2 6 


(2.13) 


As stated above V vv > should be derived from the static potential seen by the scattered 
electron, given by Eq. 2.5. We have shown 6 that the use of an MCSCF in (2.14) 
leads to a slightly less attractive potential than one obtained from an SCF approximation. 
That result is in qualitative accord with that found by Rumble et al. 8 Very recently, 
however, Meyer et al. 9 have carefully examined the II 3 resonance in the static exchange 
approximation and have found that the resonance, using an MCSCF ground state, when the 
consistent exchange terms are included, corresponds to a slightly more attractive effect 
than the corresponding consistent exchange approximation with an SCF ground state. 
Their conclusion (which we find quite surprising, but do not question) is that the use of 
SCF orbitals for exchange in Eqs. 2.10 and 2.11 is sufficiently inconsistent with the use of 
an MCSCF wave function in calculating the direct term, that it gives the opposite effect 
from what is obtained by using the MCSCF wavefunction. 9 

Nevertheless, with regard to our present calculation, we believe this effect is included 
by the way polarization has been incorporated. To our static potential, we add a polar- 
ization potential, as described in Ref. 6: 


V(r-R) V st a.tic(r-, R) + V po ,(r; R) (2.14) 

where 

W;R) = (l - V° r »(r ;R) (2.15) 

V pol is calculated from a quasi- ab initio polarized orbital derivation of the polarization 
potential 4 , but, as seen in (2.15), the latter is diminished in magnitude by a tuning factor, 

(1 - exp[(r/r 0 )] 2 ), in which r 0 was adjusted to give the 11^ resonance at the correct (i.e. 
experimental) energy. 

We found specifically in Ref. 6 that a value of r 0 = 2.430 was required when using 
Vstatic{MCSCF ), whereas r 0 = 2.934 was required for an SCF wavefunction. (There is 
an unfortunate typographical error in Ref. 6, which read r 0 = 2.394.) What this implies, 
in agreement with Meyer et al. 9 , is that although (symbolically) 

(y, tatic + V exc h) MCSCF < ( ^static 4“ V eX ch) SCF ^ (2.16) 

because 


(Vpol) 


MCSCF 


> (Vpol) SCF ’ 


(2.17) 
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it is perfectly possible that 


(K(afic + Vexch + Vpol) MCSCF ~ (Static + Kxc* + Vpol) SCF , (2-18) 

[To repeat, the above equations culminating in Eq. 2.18 are intended as a measure of “size” 
in an average sense; in detail, both the r and R dependence of each side will be noticeably 
different from the other. We expect therefore that the left hand side, which derives from 
the better (MCSCF) wave function, does represent a considerable improvement, and that 
is what is used here.] 

Thus, from the pragmatic point of view of this calculation, we consider the criticism 
of Ref. 9 to have been overcome; in fact, we believe that the results we shall present are, 
in the region of the center of the resonance (1.6 < k 2 < 3eV and v < 2), the most accurate 
that have thus far been calculated. 

Moreover, the conclusion of Ref. 9 notwithstanding, it is completely possible that 
if one had included the effect of static polarization, exchange, and correlation simultane- 
ously, (part of such effects would be described by the exchange-polarization terms in a full 
polarized orbital treatment 10 , for example), and if one had then isolated the piece labelled 
Vstatic + Vexch separately, then it might have led to the result of Weatherford et al. 6 , that 

the MCSCF result would be less attractive than the corresponding SCF result [i.e. the 
reverse of (2.16)]. 


III. RESULTS AND COMPARISONS 


To the one resonant partial wave (II ff ), we add the four most important non-resonant 
partial waves ('E,, 1 E U , J II ,,, 1 A g ), calculated in the fixed-plus adiabatic-nuclei approxi- 
mations, as discussed in Ref. 6. Let us first show and discuss the total cross section a r 
(the sum of all energetically allowed vibrational channels, summed and averaged over ro- 
tational states in the usual way). The resonance with its famous substructure 11 , is usually 
compared, as it is here, with the experimental result of Kennedy. 12 The theoretical curves 
shown are our previous 10 state result 7 (dashed curve), our present 15 state result (solid 

curve), compared to experiment and the Schwinger multichannel calculation of Huo et 
al. 13 - 14 

With reference to our calculations, the comparison of the 10 and 15 state results gives 
a good idea of the convergence of the close coupling expansions: we would say that our 
calculations are well converged to just beyond the first resonance (it 2 < 1.95eU), reasonbly 
well converged to just beyond the second (« 2.2eV), and approximately converged to 
~ 2 -5eV. Of particular note, therefore, is the fact that ours is the only calculation which 
describes the magnitude of the first peak (a T « 27A 2 at k 2 = 1.95eK) and gives the ratio of 
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the first two peaks accurately. The calculation of Huo et al. 14 is almost as satisfactory, while 
their adjoining paper 13 gives the whole sequence of vibrational excitation cross sections in 
remarkable accord with experiment. 15 

The ability of such theories as in Refs. 13 and 14 to achieve such an elaborate 
overview of the entire resonance substructure goes back to the physical ideas underlying 
the boomerang model (cf. Ref. 17), which have been given their most rigorous justification 
in the R-matrix theory and calculations of Schneider, LeDourneuf, and VoKyLan. 18 More 
discussion of these theories will be included in the latter parts of this paper. 

In this paper, we shall concentrate on angular distributions. In Fig. 2, we show the 
elastic differential cross section at 1.50 eV, just below the onset of the U g resonance. Also 
shown are various other theoretical 2 ' 14 ' 19 and experimental 20 - 22 results. Both similarities 
and differences are evident. The results continue into the heart of the resonance region 
(k 2 ss 2. leR), in Figs. 3 and 4. We have divided those results into two parts: in Fig. 
3, we show our present results at three energies surrounding 2.10 eV; one sees how sig- 
nificantly the shape varies over 0.1 eV. This is particularly relevant to the fact that the 
latest experiment 22 only claims an energy uncertainty of just that amount. Nevertheless, 
the discrepancy between all these results and the original hybrid theory calculation 2 is 
clear. In Fig. 4, the same experimental information is compared to the present calculation 
and that of Huo et al. 14 at 2.10 eV. Here the similarity of the calculated results, both of 
which provide absolute values, is the most striking feature. It should be noted that an 
SCF target representation was used in Ref. 14. The agreement between theories suggests 
that the recent experiment 22 is dominated by the particular energy in the composite beam 
which gives the dominant cross section at a particular angle. 

The comparison of experiment and theoty at 3.0 eV is shown in Fig. 5. The similarity 
of the recently calculated results continues, but — referring back to Fig. 1 — we emphasize 
that at 3 eV, the energy is definitely pressing the outer edge of reliability of the present 
calculation. The other interesting feature of the latest experiment 22 , noted there as well, is 
that the absolute magnitude in the directions 0 < 110° favors the original hybrid result. 2 

In the next three figures (6,7,8), we show differential cross sections in the excitation 
transition to the first vibrational state. At the lowest energy (1.50 eV), our present calcu- 
lation (Fig. 6) is definitely favored by the experiment of Brennan et al. 22 , whereas, in the 
vicinity of the dominant peak (Fig. 7), that is only true in the middle of the angular range. 
Note that a different set of experimental results from an Australian group has also been 
included. 23 At the highest energy, 3 eV, the graph also includes results of an R-matrix 
calculation 24 and yet another experiment. 25 Here the similarity of the present results with 
those of the R-matrix calculation is the most notable feature, particulary at middle angles, 
where they agree best with the experimental results of Ref. 23. 

Finally we show in Fig. 9, the angular distribution associated with the excitation of the 
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second vibrational state in the vicinity of the dominant peak. Here the two experimental 
results 22 - 23 are in agreement with each other to within their experimental errors, and they 
agree best with our present result at the experimental energy 2.10 eV, rather than at its 
fringes (2.05, 2.15 eV). 


IV. DISCUSSION 


This completes the presentation of the results of this calculation. (More results are 
available upon request.) Although the method (the hybrid theory) is ab initio in principle, 
it contains here one phenomenological parameter, the polarization cut-off r 0 [cf. Eq. 2.15], 
The main practical purpose for doing this is to provide the best cross sections for (several 
space) applications. It will be recalled that our (hybrid) theory was developed in the 
context of SAR arcs. 2 The numerical results of that calculation were collected as a NASA 
document. 26 The present (more accurate) cross sections are intended for the understanding 
of secondary electron flux in the lower F region of the ionosphere. The specific question 
concerns whether or not there is a dip in the electron distribution function. 27 Relevant 
ionospheric calculations are now ongoing at Goddard Space Flight Center 28 , using our 
cross sections and those of others, principally Huo et al. 13,14 

From a more fundamental point of view, the present results are intended to be a more 
definitive comparison with experiment in the resonance region, particularly with angular 
distributions. As we have seen from the comparisons in the previous section, there is still 
insufficient agreement among experiments themselves to provide a definitive check at this 
time. In addition, future experiments will require an even finer energy resolution, at a 
finer energy grid, to be compelling in this regard. This having been said, it is important 
to acknowledge that great progress in experimental angular distributions, culminating in 
the recent work of the Australian group 22,23 , has already been made. 

With regard to calculational methodology, specifically hybrid vs. R- matrix and/or 
other L basis set theories, in addition to what was said above and elsewhere 29 , it is clear 
that the latter are capable of giving the greater overall accuracy at the present time, as is 
exemplified by the results of Refs. (13,14,18). (In further detail, for example, the R-matrix 
theory shows that the U g resonance is dominated by short range correlations rather than 
the long range polarizability.) However, if insufficient correlation is included, the R-matrix 
method and L 2 methods in general, can yield noticeable inaccuracy (cf. Ref. 30, for 
example). 

Furthermore, the basic tenet of the hybrid theory is common to both approaches; 2 ’ 29 
it is the fact that if the interaction time of the resonance is comparable to vibrational time 
scales of the target molecule, but short compared to rotation time scales, then in one way 
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or another, a dynamical treatment of the iV + l electron system is required. That treatment 
can be either a dynamical coupling of the incident particle with the vibrational motion of 
the target, or a recalculation of the Born-Oppenheimer problem of the (N -f- l)-electron 
system, followed by a calculation of the vibrational spectrum in the (necessarily complex) 
potential energy well of the compound system. 

The drawback of the present methodology is due to the slow convergence of the vi- 
brational close coupling expansion. In principle, that can be overcome by going to a 3d 
approach in which the internuclear separation ( R ) becomes the third dynamical variable. 
Such a theory has already been outlined, 31 and the non-iterative pde method 1 has been 
generalized to three and higher dimensional equations and applied to a solvable model. 32 
The implementation of this program is already in progress, however, its completion will 
not be easy: in particular, the detailed substructure that such calculations reveal, will 
depend critically on the range and mesh size with which the R variable can be covered. 

Perhaps the most exciting potential application of the 3d pde technique, also dis- 
cussed in Ref. 32, is the fact that it can be applied to scattering (in principle) from 
arbitrary polyatomics in the fixed-nuclei approximation. 33 Augmented by the adiabatic- 
nuclei approximation, 33 such a method would be an invaluable tool in studying scattering 
processes in galactic environments, such as the Orion nebula, where it is known that exotic 
molecules can form (cf. Ref. 34, for example) which may not be amenable to laboratory 
experimentation, so that theoretical calculation provides the only reasonable alternative. 
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Figure Captions 


Fig. 1. Total integrated cross section in square Angstroms (A 2 ): Experiment: (* * *) Ref. 

12. Theory: ( ) present [15 state results]; (---) our 10 state results (Ref. 7); 

(-o-o-o-) Ref. 14. 

Fig. 2. dooo/df) at 1.50 eV: Theory: ( ) present; (-o-o-o-) Ref. 14; ( ) Ref. 2; 

( ) Ref. 19. Experiment: ($$$) Ref. 22; (+++) Ref. 20; (AAA) Ref. 21. 

Fig. 3. dcroo/dfi. Theory: ( ) present at 2.05 eV; ( ) present at 2.10 eV; 

( ) present at 2.15 eV; ( ) Ref. 2. Experiment: ($$$) R e f. 22; (+ + + ) 

Ref. 20. 


Fig. 4. 

dcroo/dS7 a t 2.10 eV. Theory: ( 

— ) present; (-o-o- 

-o-) Ref. 14. 

Experiment: 


($$$) Ref. 22; (+++) Ref. 20. 




Fig. 5. 

d<7 00 /dft at 3.00 eV. Theory: ( — 

) present; (- - 

-) Ref. 2. 

Experiment: 


($$$) Ref. 22; (+++) Ref. 20. 




Fig. 6. 

dooi/dfl at 1.50 eV. Theory: ( — 

) present; (- - 

-) Ref. 2. 

Experiment: 


($$$) Ref. 22; (AAA) Ref. 21. 

Fig. 7. da 0 i/dft: Theory: ( ) present at 2.05 eV; ( ) present at 2.10 eV; 

( ) present at 2.15 eV; ( ) Ref. 2. Experiment: ($$$) Ref. 22; (ooo) 

Ref. 23. 

Fig. 8. d<7oi/dfl at 3.00 eV. Theory: ( ) present; ( ) Ref. 2; ( ) Ref. 24. 

Experiment: ($<M) Ref. 22; (ooo) Ref. 23; (•••) Ref. 25. 

Fig. 9. da 02 /dV . : Theory: ( ) present at 2.05 eV; ( ) present at 2.10 eV; 

( ) present at 2.15 eV; ( ) Ref. 2. Experiment: (<f>4><f>) Ref. 22; (ooo) 

Ref. 23. 





oh- cr_j cjco — 


ENERBY(eU) 


CROSS SECTION (Ang**2) 


e“+N2 : v(0~>0) 

DIFFERENTIAL C/S 
1.50 eV 



e'+N 2 : v(0~>0) 

DIFFERENTIAL C/S 
2.05, 2.10, 2.15 eV 




e"+N2: v(0~>0) 

DIFFERENTIAL C/S 
2.10 eV 



ANGLE (deg) 



e"+N2 : v(0— >0) 

DIFFERENTIAL C/S 
3.00 eV 



CROSS SECTION (10**(-18)cm**2/sr) 


e"+N2 : v(0~>l) 

DIFFERENTIAL C/S 
1.50 eV 



ANGLE (deg) 


CROSS SECTION (10**(-18) cm**2/sr) 


e'+N 2 : v(0— >1) 

DIFFERENTIAL C/S 
2.05, 2.10, 2.15 eV 



ANGLE (deg) 



e"+N2 : v(0-->l) 

DIFFERENTIAL C/S 
3.00 eV 



ANGLE (deg) 


e"+N2 : v(0— >2) 

DIFFERENTIAL C/S 
2.05, 2.10, 2.15 eV 



ANGLE (deg) 


